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ABSTRACT 


An  efficient  method  of  solving  the  three-body  Schroedinger  equation  Is 
presented.  The  wave  function  is  decomposed  into  the  product  of  a 
correlation  factor  describing  the  singularity  and  clustering  structure,  and 
a  smooth  factor  expanded  In  hyperspherlcal  harmonics.  The  application  to 
the  Helium  atom  yields  a  ground  state  energy  of  2.9037244  (2.5033052)  au  for 

Infinite  (finite)  nuclear  mass.  The  convergence  pattern  shows  that  the 
accuracy  of  these  values  Is  better  than  a  few  parts  In 


The  hypersphen'cal  harmonic  method1"5  is  one  of  the  most  general 
methods  of  solving  the  few-body  Schroedinger  equation.  Theoretically,  this 
method  is  able  to  handle  problems  involving  three,  four,  five  or  more 
interacting  particles.  Practically,  however,  the  huge  degeneracy  of  the 
hypersphen'cal  basis  and  the  slow  convergence  of  the  method  for  realistic 
interactions  lead  to  intimidatingly  large  sets  of  coupled  (one-dimensional) 


equations  even  for  modern  computers  to  handle,  u.erefore,  most  calculations 
to  date  have  either  employed  only  a  limited  number  of  hypersphericals5"1^, 
which  leads  to  a  significant  loss  in  accuracy,  or  have  employed  a  select. 


restricted  basis,  e,g.,  the  "optimal  subset  method"  used  by  Fabre  and 
others.11"15 


Clearly,  when  the  wave  function  is  expanded  in  a  complete  set  of 
hypersphen'cal  basis  functions,  one  is  assured  that  the  wave  function  and 
expectation  values  converge  to  the  correct  values.  This  advantage  is  not 
present  in  ad  hoc  variational  wave  functions  or  In  wave  functions  built  on 
an  incomplete  set.  The  main  purpose  of  the  present  paper  is  to  introduce  a 
technique  to  dramatically  hasten  the  convergence  of  the  hyperspherical 
expansion  in  order  to  make  it  suitable  for  high  precision  calculations. 


where  X  is  known  and  chosen  to  take  into  account  singularities  and 
clustering  properties  of  the  wave  function,  and <p  is  the  part  to  be  expanded 
in  hyperspherical s.  Of  course,  any  wave  function  can  be  expressed  in  the 
form  of  Eq.  (1),  so  for  the  method  to  be  practical  X  roust  be  reasonably 
simple  to  obtain,  the  equations  for  ^  should  not  be  overly  complicated, 
and  tp  should  be  smooth  enough  so  as  to  be  readily  expressible  in  a  small 
number  of  hyperspherical s.  We  will  demonstrate  the  fulfillment  of  these 
conditions  for  the  case  of  the  ground  state  of  the  helium  atom. 

The  idea  to  separate  a  "configuration"  part  X  of  the  wave  function  to 

improve  the  convergence  of  the  hyperspherical  harmonic  expansion  is  not 

new.  A  separation  of  the  form  has  been  employed  to  apply  the 

16  1  fi 

hyperspherical  technique  in  nuclear  scattering  problems  ,  but  these 

calculations  either  only  employ  a  very  few  hyperspherical s  or  do  not  use 

realistic  nuclear  forces.  A  separation  of  the  product  form  of  Eq.  (1)  has 

19 

been  applied  to  nuclear  bound  states  ,  but  here  the  optimal  subset  and 
other  approximations  were  employed.  Furthermore,  predicted  nuclear 
properties  were  compared  to  those  calculated  by  completely  different  other 
methods,  which  were  themselves  approximate.  Thus  any  firm  assessment  of  the 
ultimate  accuracy  of  the  decomposition  of  Eq.  (1)  was  difficult  to  make. 

To  obtain  a  better  understanding  of  the  product  decomposition  (1),  we 

limit  ourselves  in  this  paper  to  the  problem  of  three  Coulomb-bound 

particles  -  specificrMy  the  Helium  atom.  Here  the  potential  is  exactly 

known  and  many  very  precise  calculations  have  been  carried  out  by 
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variational  methods,  so  comparisons  can  be  done.  More  importantly, 
the  Coulomb  potential  has  exactly  the  same  type  of  singularity  at  the  origin 
as  realistic  nuclear  potentials  (i.e.,  Yukawa  e  _/<r/r),  and  it  is  the 


nature  of  this  singularity  that  leads  to  the  slow  convergence  of  the 
hyperspherlcal  expansion  In  either  atomic  or  nuclear  systems.  Indeed,  for 
the  1/r  type  of  behavior,  the  hyperspherlcal  expansion  for  large  maximum 
global  angular  momentum  Km  converges  as  a/K^  for  wave  functions 
and  b/K^  for  energies,  compared  to  the  exponential  convergence 
(*/e‘cKm)  for  smooth  (e.g.,  Gaussian)  potentials10’13,14,24,25.  While 
the  power  or  exponential  laws  of  convergence  depend  on  the  singularity 
structure  of  the  potential,  the  coefficients  a,b,c.  Involved  in  the 
convergence  laws  ,  which  also  determine  the  accuracy  of  the  expansion, 
depend  on  other  features  of  the  system,  most  notably  the  cluster  structure. 
In  the  case  of  clustered  nonsymmetric  dispostion  of  particles  In  space, 
these  coefficients  tend  to  be  large  leading  to  significant  discrepancies  of 

c 

results  for  finite  K  with  "exact"  results.  Indeed,  it  Is  known  that 

m 

one  must  include  Km=KR  global  angular  momentum  states  to  accurately 
describe  a  system  of  binding  energy  K2/2M  (M  3  particle  mass)  and 
characteristic  dimension  R.  Hence,  in  the  case  of  an  excited  atom  in  which 
a  clustered  configuration  exists  where  one  (or  more)  constituents  is 
positioned  far  away  from  the  rest,  convergence  will  be  slow  even  for  smooth, 
nonsingular  potentials. 

Since  convergence  Is  determined  by  both  singularity  and  cluster 

properties,  the  hope  is  that  a  simple  but  clever  choice  ofX  ,  reflecting 

the  physical  peculiarities  of  the  problem,  will  yield  a  smooth  function 

describable  with  only  a  few  hyperspherlcal  harmonics.  One  simple  choice 
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for  X  is  that  of  product  correlation  functions  of  the  Jastrow  type 


where  r. .  Is  the  distance  between  particles  1,  j.  Since  for  the  bound 
1 J  .  ’ 

state  the’wave  function  in  (1)  should  be  exponentially  decreasing  with 
interparticle  distance,  a  natural  choice  of  f  could  be 


£ 


■g  Ly 


e  "‘S 


(3) 


where  Tf. .  are  adjustable  parameters.  We  shall  shortly  see  that  this 

*  J 

exponential  type  is  also  appropriate  to  handle  the  short-range  correlations 
brought  about  by  1/r^  singularities  in  the  potential. 


In  the  case  of  the  ground  state  of  the  helium  atom,  the  spatial  wave 
function  must  be  symmetric  under  interchange  of  the  two  electrons.  It  is 
convenient  to  choose  X  and  separately  symmetric,  which  means 

‘  S0  ^  h3Ve 

(4) 

where  £=  3^.  Substitution  of  (1)  and  (4)  into  the  three-body 
Schroedinger  equation,  in  the  center-of-mass  frame 


where 

yields 


(5a) 

(5b) 


(6a) 

M'<pd ,%)=  (-±vx+v')p(i?;):E$(Zv, l6b> 


where 


and  whore  the  effective  (velocity-dependent)  potential  V  is  given  by 


V':  (‘VA*X>V-jL'j'Xtf-2/ra-2/ra3+i/i;  .  <7) 


me  =  11  =  1 , 


Eqs.  (5)  -  (7)  are  written  in  atomic  units  (au),  where  e 
with  Jacobi  coordinates  -  r2),|=/ 9  (r1  +  rg  -2r3), 

=  1/  ,/T, />  =  c*  Jm/M+2 ,  with  Z,  M  being  the  nuclear  charge  and  mass, 
respectively.  The  operator  ^  is  the  six-dimensional  gradient  operator 
given  by  V  =  ). 


One  choice  of  3*  and  <5  that  reduces  the  degree  of  singularity  of  the 
Coulomb  potential  (or  any  l/r.^  potential)  at  zero  interparticle  distance 
is  the  choice 


2T=  M  2  /(Air/)  ,  S--I/Z.  '8> 


In  this  case,  as  one  can  easily  check,  the  1/r  type  singularity  In  *  ie 

Coulomb  potential  are  exactly  cancelled  by  corresponding  singularities  in 

the  -  1/2  XJX)C/X  term,  so  the  effective  potential  is  never 

✓ 

infinite.  However,  V  contains  scalar  products  of  unit  vectors  along 
interparticle  distances.  These  scalar  products  suffer  finite 
discontinuities  when  any  zero  interparticle  distance  (r..  =  0)  Is  crossed 

*  J 

in  any  direction,  so  the  effective  potential,  while  discontinuous  at 
r^j  *  0,  at  least  Is  bounded.  The  result  Is  that  the  second  derivative 
of  ^  Is  bounded  (though  discontinuous)  with  an  Infinite  third  derivative. 
In  the  usual  Schroedlnger  equation  with  the  Coulomb  potential  it  Is  the 
first  derivative 


that  is  bounded  and  discontinuous  (i.e.,  there  are  cusps  in  the  wave 

function).  This  means  the  hyperspherical  expansion  for  ^  will  converge 

3  2 

as  ~  1/Kq  as  opposed  to  ^  1/K^. 

Note  that  the  -  ( )'7  term  in  the  effective  potential  in  (7) 
is  nonhermitian.  Equations  (5)  and  (6)  are,  nevertheless,  strfctly 
equivalent  and  have  identical  eigenvalues.  Approximation  of  function^*  by  a 
finite  sum  of  hyperspherical  harmonics,  however,  destroys  the  equivalence. 
The  eigenvalues  of  the  truncated  versions  of  (6)  could,  in  principle,  even 
be  complex.  However,  they  approach  real  values  both  in  the  limits  of  just 
one  hyperspherical  harmonic  and  when  the  set  becomes  complete.  Since 
eigenvalues  are  real  in  both  of  these  disparate  limits  {N  =  1,  M +oo  ),  one 
should  not  be  surprised  that,  in  practice,  the  eigenvalues  of  equation  (6) 
remain  real  for  truncation  at  any  number  of  hyperspherical  functions.  This 
is,  in  fact,  the  case  in  our  actual  calculations. 

Another  consequence  of  the  nonhermiticity  of  the  Hamiltonian  of  Eq.  (6) 
is  that  the  error  in  the  eigenvalue  E  of  this  equation  is  no  longer 
proportional  to  the  square  of  the  error  in  the  wave  function.  This  property 
is  applicable  to  the  usual  hyperspherical  expansion  of  the  wavefunction  of 
Eq.  (5),  but  not  to  our  modification  which  instead  expands  solutions  of  Eq. 
(6).  Therefore,  we  cannot  expect  the  energy  eigenvalue  of  Eq.  (6)  to 
necessarily  converge  monotonically  from  above.  It  now  may  oscillate  around 
the  correct  limit.  More  precisely,  since  for  a  normalized  wavefunction  ^ 


we  should  expect,  if  E  is  real,  that  6e=E-E^  =  . 

Since  Vk£  ,  while  S^X1)"  1/4  due  to  the  fact  that 
2 

contains  cusps,  which  reflects  the  singularities  of  the  Coulomb 
potential ,  £ E  converges  as  1/K|j|.  That,  means  that  differences  3 

,«V 

-  Ek  2  between  energy  values  calculated  for  different  K  fall 
off  at  least  as  fast  as  1/K^,  although  nonmonotonically. 

The  monotonic  convergence  of  1/K^  can  be  recovered  by 

equating  the  energy  to  the  expectation  value  of  the  true  Hamiltonian 
calculated  with  the  wave  functions  If'  obtained  in  our  method  instead  of  to 
the  eigenvalues  of  Eg.  (6).  (The  two  prescriptions  are  not  equivalent  as 
they  a»*e  in  the  usual  hyperspherical  treatment).  In  this  letter  we  present 
results  for  the  eigenvalues  as  we  have  not  as  of  yet  implemented  the 
evaluation  of  wave  functions  and  expectation  values  in  our  calculations. 

The  results  of  the  calculations  of  the  ground-state  eigenvalue  of 

fq.  (<j)  for  the  helium  atom  are  summarized  in  Table  I,  along  with 

comparisons  with  the  usual  hyperspherical  method  and  results  obtained 
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by  variational  methods  .  The  numerical  method  Is  essentially  the  same 
as  in  refs.  9-10,  except  that  different  Taylor  series  are  employed  on 
numerous  sections  of  the  arid  in  the  hyperradius as  opposed  to  a  single 
MacLaurin  series  around  p  ~  0  in  these  earlier  works.  This  procedure 
enables  us  to  obtain  twelve  decimal  place  accuracy  in  the  integration  of  the 
resulting  coupled  differential  equations  with  CRAY  single-precision 
arithmetic.  Analytic  expressions  are  available  for  the  matrix  elements 
of  the  Coulomb  part  of  the  effective  potential  (7).  We  have  derived 
analytic  expressions  for 


nn$t  of  the  matrix  elements  of  the  remaining  portion  of  (7),  but  a  part  of 
the  term  quires  numerical  evaluation.  Here  using  a  16C0  by  1600 

mesh  in  the  A  and  A  variables  leads  to  an  trror  of  less  than  one 

C) 

part  in  2.9  x  10  for  the  eigenvalue  of  Eq.  (6). 

The  energy  differences  =  E^  -  EK  _2  are  given  In  Table 

II.  In  order  to  find  the  law  of  convergence,  we  have  fit  £ek  (for 

Kr  ^16)  with  exponetially  and  inverse  power  decaying  oscillating 

functions;  however,  the  best  fits  consistently  preferred  the  inverse  power 

o 

over  the  exponential  decay.  The  optimum  form,  as  judged  by  ^  per  degree 
of  freedom,  is  given  by 

S.A'  cos 

<* 

with  H  =  2,  (^  =  0  giving  the  best  fit  .  When  all  remaining  parameters, 
including  the  overall  powers  of  K~\  were  allowed  to  vary  freely,  we 
obtained  the  following  values:  AQ  =-0.04,  Aj  =  -1.98,  A2  =  0.49, = 

0.78.  The  calculated  in  this  fit  also  appear  in  Table  II.  The 

value  for  a(  ,  including  standard  error,  is  o<=  5.96  +,  .22.  This  agrees 
fairly  well  with  the  theoretically  predicted  1/K^  convergence. 

In  summary,  we  have  presented  a  method  of  direct  (nonvariational ) 
solution  of  the  three-body  Schroedinger  equation  based  on  the  expansion  of 
the  smooth  part  of  the  product  wave  function  into  a  rapidly  convergent 
hyperspherical  series.  The  application  of  the  method  to  the  calculation  of 
the  ground  state  energies  of  the  Helium  atom  shows  its  competiveness  with 
the  best  available  variational  calculations6®"6"*,  and  its  superiority  to 
other  attempts  of  direct  solutions  of  the  Schroedinger  equation,  such  as 


-  K 

K  rr 


those  based  on  the  finite  difference^  and  the  finite  element^ 
nethods.  The  possibility  of  accurate  calculations  without  assumptions  about 
the  form  of  the  wavefunction,  as  well  as  the  nonpert.urbati ve  treatment  of 
the  center-of-mass  motion,  makes  our  direct  method  of  solution  especially 
attractive  for  excited  systems  and  systems  of  particles  of  comparable 
masses.  These  systems  are  difficult  to  treat  by  the  variational  method. 

One  of  the  authors  (V.B.M.)  wishes  to  express  his  gratitude  to  the  U.S. 
pppartnent  of  Energy  for  their  support  of  this  work.  ;.'e  also  thank  Drs.  E.F. 
Pedish  and  John  Morgan  III  for  helpful  remarks. 
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Table  i 


Calculated  ground  state  energies  of  the  Helium  atom  in  au.  is  the 
maximum  global  angular  momentum  and  N  the  number  of  resulting  equations 
solved.  The  energy  values  in  the  three  columns  are  for  an  infinite  mass 
nucleus,  a  finite  mass  nucleus  ( Me/Mnucl eus  =  T*3709337  x  10’4),  and 

o  IQ 

for  the  usual  hyperspherical  expansion  *  (infinite  mass  nucleus). 


Km 

N 

(Infinite  mass) 

EK*n 

(Finite  mass) 

EIU 

(Unmodified  hyperspherical s ) 

0 

1 

2 .9767604 

2.9765759 

2.50002 

4 

4 

2.91C26686 

2.90989022 

2.78437 

8 

9 

2.90321786 

2.90281192 

2.85022 

12 

16 

2.90390765 

2.90349598 

2.87601 

16 

25 

2.90370440 

2.90328900 

2.88754 

20 

36 

2.90374150 

2.90332453 

2.89358 

24 

49 

2.90372303 

2.90330507 

28 

64 

2.90372743 

2.9033C892 

32 

81 

2.903724364 

2.903305480 

36 

100 

2.903725184 

2.903306065 

40 

121 

2.903724448 

2.903305161 

Other 

2.903723721 

2.903303721 

Methods 

2.90372437620 

2.90330436520 

2.90372437722*23 

2.90330437420 

Table  II 


The  ground  state  energy  differences  =  E^  -E^  _0  for  an 

.*^V  *■%  ‘V’V 

infinite  mass  nucleus.  £  E^  (fit)  are  the  energy  difference  obtained  in 
our  fitting  analysis  with  the  parameters  given  in  the  text. 


Kn 

£fk  (au) 

Sek  (fit)  (au) 

**v 

2 

-.093047? 

4 

.02691113 

6 

-.00375277 

8 

-.00329624 

10 

.00009660 

12 

.00059320 

14 

-.00009451 

16 

-.00010874 

-.00010205 

18 

-.00000258 

-.00001414 

20 

.00003968 

.00004279 

22 

-.00000839 

-.00000629 

24 

-.00001008 

-.00000911 

26 

-.00000216 

-.00000142 

28 

.00000656 

.00000575 

30 

-.00000145 

-.00000105 

3? 

-.000001622 

-.000001639 

34 

-.000000375 

-.000000255 

36 

.000001195 

.000001282 

38 

-.000000370 

-.000000271 

40 

-.000000366 

-.000000434 

